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Abstract 



1 \ In principle one can test the validity of charge symmetry for parton distri- 

r~| ! butions by comparing structure functions measured in neutrino and charged 

lepton deep inelastic scattering. New experiments make such tests possible; 

X 

they provide rather tight upper limits on parton charge symmetry violation 
[CSV] for intermediate Bjorken x, but appear to show evidence for CSV effects 
at small x. We examine two effects which might account for this experimental 
discrepancy: nuclear shadowing corrections for neutrinos, and strange quark 
contributions s(x) ^ s(x). We show that neither of these two corrections re- 
moves the experimental discrepancy between the structure functions. We are 
therefore forced to consider the possibility of a surprisingly large CSV effect 
in the nucleon sea quark distributions. 
PACS: 13.60.Hb, 13.15. +g, 12.40.Vv, 11.40.Ha 
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In nuclear physics charge symmetry, which interchanges protons and neutrons (simul- 
taneously interchanging up and down quarks), is respected to a high degree of precision. 
Most low-energy tests of charge symmetry find that it is good to at least 1% in reaction 
amplitudes 0. Therefore, charge symmetry is usually assumed to be valid in discussions of 
strong interactions. Currently all phenomenological analyses describe deep inelastic scatter- 
ing [DIS] data using charge symmetric parton distributions. Until recently this assumption 
seemed to be justified, since the quantitative evidence which could be extracted from high 
energy experiments, although not particularly precise, was consistent with charge symmetric 
parton distributions |J. 

Experimental verification of charge symmetry is difficult, partly because the relative 
charge symmetry violation (CSV) effects are expected to be small, requiring high precision 
experiments, and partly because CSV often mixes with parton flavor symmetry violation 
(FSV). Recent experimental measurements by the NMC Collaboration ||, demonstrating 
the violation of the Gottfried sum rule [[|, have been widely interpreted as evidence for 
what is termed SU(2) FSV. The measurement of the ratio of Drell-Yan cross sections in 
proton-deuteron and proton-proton scattering, first by the NA51-Collaboration at CERN 
and more recently by the E866 experiment at FNAL ||, also indicate substantial FSV. 
However, as pointed out by Ma M , both of these experiments could in principle be explained 
by sufficiently large CSV effects, even in the limit of exact flavor symmetry. In view of these 
ambiguities in the interpretation of current experimental data, it would be highly desirable 
to have experiments which separate CSV from FSV. A few experiments have been already 
proposed || and could be carried out in the near future. 

At the level of parton distributions charge symmetry implies the equivalence between up 
(down) quark distributions in the proton and down (up) quarks in the neutron. We define 
charge symmetry violating distributions 



where the superscripts p and n refer to quark distributions in the proton and neutron, 
respectively (quark distributions without subscripts will refer to the proton). The relations 
for CSV in antiquark distributions are analogous. Exact charge symmetry would require 
that the quantities Su(x) and Sd(x) vanish. 

In the quark-parton model the structure functions measured in neutrino, antineutrino and 
charged lepton DIS on an iso-scalar target, Nq, are given in terms of the parton distribution 
functions and the charge symmetry violating distributions || 



Su(x) = u p (x) - d n (x) 
5d{x) = d p {x) - u n {x) 



(1) 



F^(x,Q 2 ) 
F 2 m °(x,Q 2 ) 
xF» N °(x,Q 2 ) 
xF* N °(x,Q 2 ) 



x[u(x) + u(x) + d(x) + d(x) + 2s(x) + 2c(x) — 8u(x) — 5d(x)] 
x[u(x) + u(x) + d(x) + d(x) + 2s(x) + 2c(x) — 5d(x) — 8u{x)\ 
x[u(x) + d(x) — u(x) — d(x) + 2s(x) — 2c(x) — Su(x) + 6d(x)] 
x[u(x) + d(x) — u(x) — d(x) — 2s(x) + 2c(x) — Sd(x) + Su(x)] 



F^(x,Q 2 ) 



—x[u(x) + u(x) + d(x) + d(x) + -(s(x) + s(x)) + -(c(x) + c(x)) 
18 5 5 

4 1 

-(5d(x) + 5d(x)) (5u(x) + Su(x))] 

5 5 



(2) 
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By comparing structure functions measured in neutrino and charged lepton deep inelastic 
scattering [DIS] on isoscalar targets, it is possible to test parton charge symmetry. An 
example is the "charge ratio" , which relates the neutrino structure function to the structure 
function measured in charged lepton DIS 



F? N °(x) 



1^2 (i) - x(s(x) + s(x))/6 

s(x) — s(x) 4Su(x) — 8u(x) — 45d(x) + Sd(x) 

Here, we defined Q(x) = Y, q = u d s(q( x ) + Q( x )) ~ 3(s(sc) + s(x))/B, and we have expanded 
Eq. |] to lowest order in small quantities. A deviation R c (x) ^ 1, at any value of x, must 
arise either from CSV effects or from an inequality between strange and anti-strange quark 
distributions. 

Recent experimental measurements make it possible to carry out a precise comparison be- 
tween F%{x, Q 2 ) and F%{x, Q 2 )- The CCFR-Collaboration compared the neutrino structure 
function F%(x, Q 2 ) extracted from their data on an iron target M with F^{x, Q 2 ) measured 



for the deuteron by the NMC Collaboration [10|. In the region of intermediate values of 



Bjorken x (0.1 < x < 0.4), the two structure functions are in very good agreement; in this 
region we can set upper limits on parton CSV contributions of a few percent. In the small 
x-region however (x < 0.1), the CCFR group found that the two structure functions differ 
by as much as 10-15%. This can be seen in Fig.|l| where the "charge ratio" has been obtained 
by integrating over the region of overlap in Q 2 of the two experiments. The various data 
points in Fig. [I] represent different ways of calculating nuclear shadowing corrections, as we 
will discuss. Several corrections must be applied to the data before any conclusions may 
be drawn from this discrepancy. The CCFR Collaboration made a careful study of overall 
normalization, charm threshold and iso-scalar correction effects [[|. Here we discuss the 
most important remaining effects, heavy target corrections for the neutrinos and differences 
between strange and anti strange quark distributions. 

Nuclear shadowing corrections for neutrinos are generally accounted for by using correc- 
tion factors obtained from charged lepton reactions at the same kinematic values. A priori, 
there is no reason to assume that neutrino and charged lepton heavy target corrections 
should be identical. We re-examined heavy target corrections to deep-inelastic neutrino 
scattering, focusing on the differences between neutrino and charge lepton scattering and on 
effects due to the Q 2 -dependence of shadowing for moderately large Q 2 . This work will be 
published elsewhere We used a two phase model which has been successfully applied 



to the description of shadowing in charged lepton DIS [Tj|]. In this approach, vector meson 
dominance was used to describe the low Q 2 virtual photon or W interactions, and Pomeron 
exchange was used for the approximate scaling region. In generalizing this approach to weak 
currents, the major new features are that the axial- vector current is only partially conserved, 
and that the weak current couples to axial vector mesons in addition to vector mesons. 

Using this two phase model, we calculated the shadowing corrections to the CCFR 
neutrino data and used these corrections in calculating the charge ratio R c of Eq. |3|. The 
result is shown in Fig.|I|. The solid triangles show the charge ratio when no shadowing 
corrections are used. The open circles show the charge ratio when the neutrino data is 
modified using heavy target shadowing corrections from charged lepton reactions, and the 
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solid circles show the results calculated specifically for neutrinos using our two phase model. 
For x > 0.1, the two shadowing corrections give essentially identical results. At small x, 
careful consideration of neutrino shadowing corrections decreases, but does not resolve, the 
low-x discrepancy between the CCFR and NMC data. 

At this point it is instructive to review how the structure functions are extracted in 
neutrino reactions. Because of the extremely small cross sections, it is necessary to inte- 
grate cross sections over all energies in order to accumulate sufficient flux. When this has 
been done, the resulting neutrino and antineutrino cross sections give two equations in the 
neutrino structure functions. If we assume that the neutrino and anti-neutrino structure 
functions are equal F 2 (x,Q 2 ) = F 2 (x,Q 2 ), with an analogous relation for xF 3 (x, Q 2 ), then 
we have two linear equations in two unknowns. From Eq. ^| we see that these relations will 
be true if charge symmetry is valid and if s(x) = s(x); there is an additional correction to 
xF% and xF% from strange quarks. 

The resulting structure function F 2 CFR is a flux weighted average between neutrino 
and anti- neutrino structure functions (see Ref. 0). This fact becomes important if charge 
symmetry is violated or the strange and anti-strange quark distributions are different. If we 
define a = N V /(N V + N&), where N„ and Np are the number of neutrino and anti-neutrino 
events, respectively, F 2 FR (x, Q 2 ) is proportional to 

F 2 CCFR (x, Q 2 ) = aF%(x, Q 2 ) + (1 - a)F»(x, Q 2 ) 

= l -[F»(x,Q 2 ) + F*{x,Q 2 )} + l(2a- l)[i^(*,Q 2 ) - if(x,Q 2 )]. (4) 

This is equal to \[F 2 (x, Q 2 ) + F 2 (x, Q 2 )} if a = \ or if the two structure functions are equal 
(which implies s(x) = s(x) and the validity of charge symmetry). The CCFR-Collaboration 
collected 1, 300, 000 neutrino and 270, 000 anti-neutrino events 0; thus a ~ 0.83 so that to 
a good approximation F 2 CFR (x, Q 2 ) can be regarded as a neutrino structure function. 

The most likely explanation for the small- a: discrepancy in the charge ratio is either from 
different strange quark distributions s{x) ^ s(x) or from charge symmetry violation. 



First, we examine the role played by the strange and anti-strange quark distributions. As- 
suming charge symmetry, the strange and anti-strange quark distributions are given by a 
linear combination of the structure functions measured in neutrino and in muon DIS, 

^F 2 CCFR (x, Q 2 ) - 3F 2 NMC (x, Q 2 ) = \x [s(x) + *(*)] + ~(2e* - 1) x [s(x) - s(x)}. (5) 
o 2 o 

Under the assumption s(x) = s(x), this relation could be used to extract the strange quark 
distribution. However, as is well known, s(x) obtained in this way is inconsistent with the 
distribution extracted from independent experiments. 

The strange quark distribution can be determined directly from opposite sign dimuon 
production in deep inelastic neutrino and anti-neutrino scattering. The CCFR Collaboration 
performed a LO |14| and NLO analysis |15[ of their dimuon data using the neutrino (anti- 



neutrino) events to extract the strange (anti-strange) quark distributions. They found that 
s(x) and s(x) were equal within experimental errors in NLO |TJJ. However, since the number 
of anti-neutrino events is much smaller than that of the neutrino events, the errors of this 
analysis are inevitably large. 
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Since the dimuon experiments are carried out on an iron target, shadowing corrections 
could also modify the extracted strange quark distribution. The CCFR-Collaboration nor- 
malized the dimuon cross section to the "single muon" cross section and argued that the 
heavy target correction should cancel in the ratio. However, the charm producing part of the 
structure function F£ p (x, Q 2 ) could in principle be shadowed differently from the non-charm 
producing part i^ cp (x,Q 2 ). We tested this hypothesis by calculating the neutrino shadow- 
ing corrections for both the charm and non-charm producing part of the structure function. 
The results will be published separately ||16|| . We find that, while the relative importance 
of the Pomeron and VMD components are different in the charm producing (cp) and the 
non-charm producing (ncp) parts, there is essentially no difference in the total shadowing. 

It appears plausible that the low-x discrepancy in the charge ratio of Eq. [3] can be ac- 
counted for by allowing s(x) ^ s(x). To test this hypothesis we combined the dimuon 
production data, averaged over v and v events, with the structure functions from neutrino 
and charged lepton scattering (Eq.(^])). Defining a' = N U /(N U + N p ), where N v = 5,030, 
N p = 1, 060 (a' ~ 0.83) are respectively the v and v events from the dimuon production ex- 
periment [15|], the flux- weighted experimental distribution xs{x)^ from dimuon production 
is 

xs^{x) = ^x [s(x) + s(x)} + ^(2a' -l)x [s(x) - s(x)}. (6) 

This equation together with Eq. (^) forms a pair of linear equations which can be solved for 
s(x) and s(x). We can simultaneously test the compatibility of the various experiments. 

In Fig. § we show the results obtained for xs(x) (open circles) and xs(x) (solid circles) by 
solving the resulting linear equations, Eqs. |5] and [| The results are completely unphysical, 
since the extracted anti strange quark distribution is negative. Our analysis strongly suggests 
that requiring charge symmetry, but allowing s(x) ^ s(x), cannot resolve the discrepancy 
between F2 CFR (x,Q 2 ) and F 2 NMC (x,Q 2 ). The experimental results are incompatible, even 
if s(x) is completely unconstrained JT7| . 

As neither neutrino shadowing corrections nor allowing s(x) ^ s(x) removes the low- a; 
discrepancy between the neutrino and muon structure functions, there remain two possible 
explanations. Either one of the experimental structure functions (or the strange quark 
distribution) is incorrect at low x, or parton charge symmetry is violated in this region. 
Assuming the possibility of parton CSV, we can combine the dimuon data for the strange 
quark distribution (Eq. H) with the relation between neutrino and muon structure functions, 
Eq. |5| to obtain 

5 -F 2 CGFR (x, Q 2 ) - 3F 2 ^ C (x, Q 2 ) - xs^(x) = x{2 ° ~ 1} [s(x) - s(x)} 
o 3 

+ | [(5a - l)(Sd(x) - 8u{x)) + (4 - ha){5d(x) - 5u(x))\ 
6 

x(2a — 1) 1 - 

w -^-^ ~[ s ( x ) ~ s( x )\ + 2 x l Sd ( x ) ~ 8u{x)\. (7) 

In Eq. |7] we have used the experimental value a = a'. Since the experimental discrepancy is 
primarily in the small x-region, where sea quark distributions are much larger than valence 
quarks, CSV effects should appear predominantly in the sea quark distributions. Setting 
Sq v (x) = 5q(x) — Sq(x) ~ in this region gives the second relation in Eq.(|7|). 
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The left hand side of Eq. [7| is positive. Consequently, the smallest CSV effects will be 
obtained by setting s(x) = 0. In Fig.[3]we show the magnitude of charge symmetry violation 
needed to satisfy the experimental values in Eq. [7]. The solid circles are obtained if we set 
s(x) = 0, and the open circles result from setting s(x) = s(x). Both the structure functions 
and the dimuon data have been integrated over the overlapping kinematical regions. In 
averaging the dimuon data we used the CTEQ4L parametrization for s M/i (x) |19[], and we 
extracted the strange quark distributions according to Eq. The CSV effect required to 
account for the NMC-CCFR discrepancy is extraordinarily large. It is roughly the same size 
as the strange quark distribution at small x. This CSV term is roughly 25% of the light sea 
quark distributions for x < 0.1. The existing experimental data thus appears to require a 
very surprising, and uncomfortably large, violation of parton charge symmetry at small x. 

Clearly, CSV effects of this magnitude need further experimental verification. It is hard 
to imagine how such large CSV effects are compatible with the high precision of charge 
symmetry measured at low energies. The level of CSV required is at least two orders of 
magnitude larger than theoretical estimates of charge symmetry violation [pO| , pT|| . We will 
discuss the implications of such a large violation of charge symmetry in a subsequent paper 



|T6| . Theoretical considerations suggest that 5d(x) ~ — 5u(x) pOfl . We note that with this 
sign CSV effects also require large flavor symmetry violation. At small x, our results can be 
summarized by 

5d(x) — 5u{x) » - I "' v "' ' ~' y "' I ]r(s(x) + s(x)) 




2 

5d(x) + 5u(x) « (x<0.1) (8) 

This suggests that for x < 0.1, d n (x) « 1.25 u p (x) and u n (x) w 0.75 d?{x). If CSV effects of 
this magnitude are really present, then one must include charge symmetry violating quark 
distributions in phenomeno logical models from the outset, and re-analyze the extraction of 
all parton distributions. 

This work is supported in part by the Australian Research Council, and by the National 
Science Foundation under research contract [nsf-phy9722706| . One of the authors [JTL] 



wishes to thank the Special Research Centre for the Subatomic Structure of Matter for its 
hospitality during the period when this work was carried out. 



6 



REFERENCES 



[1] G. A. Miller, B. M. K. Nefkens and I. Slaus, Phys. Rep. 194 (1990) 1; E. M. Henley 
and G. A. Miller in Mesons in Nuclei, eds M. Rho and D. H. Wilkinson (North-Holland, 
Amsterdam 1979). 

[2] J.T. Londergan and A.W. Thomas, to be published in Progress in Particle and Nuclear 
Physics, 1998. 

[3] NMC-Collaboration, P. Amaudruz et al, Phys. Rev. Lett. 66 (1991) 2712; Phys. Lett. 

B295 (1992) 159. 
[4] K. Gottfried, Phys. Rev. Lett. 18 (1967) 1174. 
[5] NA51-Collaboration, A. Baldit et al, Phys. Lett. B332 (1994) 244. 
[6] E866-Collaboration, E. A. Hawker et al, to be published in Phys. Rev. Lett. 
[7] B.-Q. Ma, Phys. Lett. B274 (1992) 433; B.-Q. Ma, A. W. Schafer and W. Greiner, 

Phys. Rev. D47 (1993) 51. 
[8] J. T. Londergan, S. M. Braendler and A. W. Thomas, Phys. Lett. B424, 185 (1998); 

J. T. Londergan, Alex Pang and A.W. Thomas, Phys. Rev D54 (1996) 3154. 
[9] CCFR-Collaboration, W.G.Seligman et al, Phys. Rev. Lett. 79, 1213 (1997) and 
W.G.Seligman Ph.D. Thesis, Nevis Report 292. 
[10] NMC-Collaboration, M.Arneodo et al, Nucl. Phys. B483, 3 (1997). 
[11] C. Boros, J. T. Londergan and A. W. Thomas, (preprint [hep-ph / 98044 1 1|) . 
[12] J. Kwiecinski and B. Badelek, Phys. Lett. B208, 508 (1988); W. Melnitchouk and A.W. 

Thomas, Phys. Lett. B317, 437 (1993); Phys. Rev. C52, 3373 (1995). 
[13] W. Melnitchouk and M. Malheiro, Phys. Rev. C55, 431 (1997); X. Ji and J. Tang, Phys. 
Lett. B362, 182 (1995); H. Holtmann, A. Szczurek and J. Speth, Nucl. Phys. A596, 
631 (1996). 

[14] S.A.Rabinowitz et al, CCFR-Collaboration, Phys. Rev. Lett. 70, 134 (1993). 
[15] CCFR-Collaboration, A.O. Bazarko et al, Z. Phys. C65, 189 (1995). 
[16] C. Boros, J.T. Londergan and A.W. Thomas, to be published. 

[17] In Ref. [|18j it was suggested that allowing s(x) ^ s(x) could account for the difference 
between the two determinations of the strange quark distribution. This result was ob- 
tained by assuming a — | in Eq. |5| and a' = 1 in Eq. ^. The choice a — | does not 
agree with experiment (a = 0.83). 

[18] S.J.Brodsky and B.Q.Ma, Phys. Lett. B381, 317 (1996). 

[19] H. L. Lai et al, Phys. Rev. D55, 1280 (1997) 

[20] E. Rodionov, A. W. Thomas and J. T. Londergan, Mod. Phys. Lett. A9, 1799 (1994); 

C.J. Benesh and J.T. Londergan, preprint nuc/-i/i/9803017. 
[21] E. Sather, Phys. Lett. B 274, 433 (1992); C.J. Benesh and T. Goldman, Phys. Rev. 

C55, 441 (1997). 



7 



FIGURES 



A 

CM 

o 



v 



1.2 



,° 1.1 



1 



0.9 



0.8 



A without shadowing 
O with shadowing (a) 
• with shadowing (b) 




Ttt 



10 



10 



FIG. 1. The "charge ratio" as a function of x calculated using the CCFR || data for neutrino 
and NMC || data for muon structure functions. The data have been integrated over the overlapping 
kinematical regions and have been corrected for heavy target effects using a parametrization (dotted 
line) for heavy target corrections extracted from charged lepton scattering (shadowing "(a)"). The 
result is shown as open circles. The ratio obtained without heavy target corrections and with 
shadowing corrections calculated in the "two phase" model (shadowing "(b)") are shown as solid 
triangles and circles, respectively. The calculated heavy target correction factors for neutrino and 
for charged lepton scattering are represented by solid and dashed lines, respectively. Statistical and 
systematic errors are added in quadrature. 
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FIG. 2. The strange quark distribution s(x) (open circles) and anti-strange distribution s(x) 
(solid circles) extracted by combining the CCFR and NMC structure functions with s(x) extracted 
from dimuon experiments, as given in Eqs. || and|(| The difference between the CCFR neutrino and 
NMC muon structure functions %F§ R ~ 3F^ is shown as solid triangles. The strange quark 
distribution extracted by the CCFR in a LO-analysis (Ref. |14|]) is shown as solid stars, that from 
a NLO-analysis (Ref. |]l5)) is represented by the solid line with a band indicating ±ler uncertainty 
in the distribution. Statistical and systematic errors are added in quadrature. 
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FIG. 3. Charge symmetry violating distributions extracted from the CCFR and NMC structure 
function data and the CCFR dimuon production data under the assumption that s(x) = s(x) (open 
circles) and s(x) « (solid circles). Statistical and systematic errors are added in quadrature. 
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